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ABSTRACT 


A  study  of  the  CABSES  vertical  plane  motions  using  the 
6-D.O.F.  Program  for  the  XR-3  craft  was  performed  with  various 
values  of  air  flow  rates  in  order  to  observe  the  sensitivity 
of  the  craft  heave,  c,  g.  acceleration  and  plenum  pressure 
to  variations  in  air  flow  rate.  The  frequency  response  due 
to  regular  and  complex  seas  as  well  as  the  transient  response 
due  to  a  weight  removal  were  used  to  observe  the  effect  of 
changes  in  air  flow  rate.  A  linear  heave-only  model  is  shown 
to  provide  useful  information  concerning  the  effect  of  air 
flow  rate  on  the  high  frequency  characteristics  of  the  XR-3 
craft . 
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I. 


INTRODUCTION 


One  important  design  consideration  of  any  naval  vessel  is 
the  habitability  of  the  craft  for  all  possible  environmental 
conditions  encountered  in  its  operation.  For  the  CAB  type  surface 
effect  ship,  it  is  important  that  the  vibrations  induced  by  the 
high  speed  operations  over  the  full  spectrum  of  the  sea  surface 
wave  profile  do  not  produce  a  great  discomfort  factor  for  the 
crew.  In  particular,  the  possibility  of  large  vertical  accelerations 
for  prolonged  periods  of  operation  should  be  avoided. 

The  primary  purpose  of  this  study  is  to  obtain  data  relating 
the  effects  of  plenum  air  flow  leakage,  and  air  influx  rates  on 
the  vertical  motions  of  The  Captured  Air  Bubble  craft.  To  do 
this  several  models  and  several  kinds  of  simulation  tests  are 
used.  The  models  used  are  a  Six  Degree  of  Freedom  computer 
model  (originally  written  by  Oceanics,  Inc.,  and  modified  for  the 
XR-3  craft  at  USNPGS) ,  and  the  simplified  linear  model  for  heave 
motions  only  (written  at  USNPGS,  see  Report  NPS62-77-002 ,  Ref.  2). 

The  simulation  tests  included  sinusoidal  sea  frequency 
response  tests,  complex  sea  frequency  response  tests,  and  step 
weight  removal  transient  tests. 
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II.  SIMPLIFIED  XR-3  MODEL 


In  order  to  investigate  the  characteristics  of  the  vertical 
motion  of  the  CAB  craft,  it  is  necessary  to  obtain  a  mathematical 
model  that  can  readily  be  used  for  analytical  purposes.  This  was 
the  principal  objective  in  Reference  1  in  which  a  simplified 
heave-only  model  of  the  XR-3  was  developed.  For  the  convenience 
of  ready  reference,  the  development  of  this  model  has  been 
repeated  in  Appendix  A. 

This  simplified  XR-3  model  has  proven  to  be  quite  useful  in 
developing  a  scaling  method  for  the  heave  motion  equations  as 
reported  in  Reference  2.  The  linear  model  used  was  described 
by  the  state  variable  vector 

x  =  (6 z,  Sz,  <SVm)T 
and  the  output  variable, 

y  =  6Pb 

as  given  in  Appendix  A  where 

Xi  =  6z,  incremental  heave  (draft) 

x 2  =  <5z,  incremental  heave  rate 

x3  =  6V  ,  incremental  air  mass  volume  with  6V  =  m,  /p 
m  m  d  a 

y  =  <$Pb,  incremental  plenum  pressure 

The  step  weight  transient  used  both  in  Reference  1  and  Reference  2 
had  the  state  equation 

x  =  A  x 
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with  initial  conditions 


x  ( 0 )  =  [6z  (0)  0  6Vm(0)]T 

which  was  brought  about  by  the  sudden  removal  of  10%  of  the  craft 
weight.  The  characteristic  equation  for  the  system 

det (si  -  A)  =  0 

is  shown  in  Appendix  A  to  be  of  third  order.  In  this  report,  the 
sinusoidal  sea  input  transfer  function  for  this  model  is  developed 
in  the  following  section.  The  step  weight  transient  is  investigated 
in  Section  IV  to  correlate  the  results  of  the  sinusoidal  sea  input 
investigation . 


III.  SINUSOIDAL  SEA  INPUT 


The  characteristics  of  the  fully  developed  sea  surface  can  be 
approximated  by  a  long  crested  wave  obtained  by  the  addition  of  a 
finite  number  of  sinusoidal  waves  of  various  amplitudes  as  described 
in  References  3  and  4  and  summarized  below. 

For  the  CAB-SES  frequency  response  analysis,  it  is  necessary  to 
develop  a  transfer  function  for  the  sinusoidal  sea  input.  What 
follows  is  a  development  of  the  transfer  function  for  the  heave  motion 
equation  of  the  simplified  XR-3  craft. 


A  Development  of  the  Simplified  XR-3  Transfer  Functions 
1.  Surface  Wave  Profile 

The  surface  wave  profile  for  both  deep  and  shallow  waves  can  be 
represented  by  a  regular  swell  or  wave: 

WAV(x,t)  =  A  sin[k  (x  +  V  t) ] 

n  w 


where:  A  =  surface  wave  amplitude  (half-height 

from  crest  to  trough)  (ft) 

k  =  the  "wave  number",  2ir/L  (ft-1) 

n  '  '  w 


L  =  wave  length 
w 


(ft) 


=  wave  velocity  or  "celerity"  (ft/sec) 


[for  deep  water,  vw  =  m  2n 


A  more  convenient  form  can  be  obtained  by  using  circular  (inci¬ 
dent)  frequency,  =  2ir/Tw  where  Tw  is  the  period,  that  is,  the 
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time  required  for  the  wave  to  travel  one  wave  length.  From  the 

Lw  /2ttL\% 

equation  for  Vw  ,  Tw  =  ^  =(-5—) 


w 


«-» «i  -  £ - 

w  \  w  / 


(k  g) 55  =  k  V 
n^  n  w 


thus  , 


WAV  ( x ,  t )  =  A  sin(knx  +  OKt) 


When  observed  at  a  fixed  point  in  space  with  x  =  0 


WAV ( 0  ,  t )  =  A  sin (w^t) 


Alternately,  if  the  wave  profile  is  observed  at  t  =  0, 


WAV(x,0)  =  A  sin(k  x) 

n 


The  slope  of  the  wave  surface  is  obtained  by  differentiation, 
WAV ( x , 


|j[wAV(x,0)]  =  kn. 


A  cos (k  x) 
n 


and  the  maximum  slope  occurs  when  knx  =  0,  it,  etc.,  thus 

.  ,  2ttA  w 

max  slope  =  k  A  =  = -  =  = - 

r  n  L  L 

w  w 

where  =  length  of  the  wave 

hw  =  2A  =  wave  height  from  hollow  to  crest. 

2.  Wave  Input  for  the  Heave  Equation 

In  the  case  of  the  CAB-SES,  the  wave  introduces  (1)  a  change  in 
the  plenum  volume  and  thus  a  change  of  plenum  pressure  and  (2)  a 
change  in  sidewall  immersed  volume  and  thus  a  change  in  the  forces 
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of  buoyancy.  The  additional  forces  due  to  the  seals  will  be  neglected 
in  this  study  since  they  largely  affect  pitch  motion  which  is  not 
considered  in  this  model.  Since  the  objective  is  to  obtain  transfer 
function  of  the  system  to  be  analyzed  for  the  condition  of  a  forcing 
function  that  is  sinusoidal  in  time,  it  is  necessary  to  recast  the 
wave  profile  equation  from  a  function  of  displacement  and  time  to  a 
function  of  time  only.  This  is  readily  possible  for  the  heave-only 
model  since  all  the  forces  along  the  hull  that  result  due  to  the  waves 
can  be  lumped  together  into  a  single  force  acting  at  the  center  of 
gravity  (CG)  which  is  the  same  point  as  the  center  of  pressure  (C.P.) 
for  this  simplified  XR-3  model. 

Given  the  wave  profile 

WAV(x,t)  =  A  sin(k  x  +  oi.t) 

n  1 

and  the  assumption  that  the  craft  is  heading  into  the  wave,  the  wave 
height,  WAV(x,t) ,  as  seen  by  the  ship  hull  is  at  the  apparent  or  en¬ 
counter  frequency,  &)e  =  ok  ( 1  +  ^— )  where  V  is  the  ship  velocity  and 

w 

is  the  velocity  of  the  wave,  thus 

WAV(x,t)  =  A  sin(k  x  +  0)  t) 

n  e 

=  A  sin  (  '°ix  +  a)  t\ 

“v -  e  ) 

w  ' 

Use  of  the  identity  for  the  sine  of  the-sum-of-two-angles  yields 

[0)  .  x  .  x 

sin(o)  t)  cos  +  cos  (co  t)sin 

e  v  e  v 

w  w 

The  differential  volume  of  the  plenum  due  to  the  wave  input  is 
d[6V  (x,t) ]  =  (W  ) [WAV (x,t) ] dx 
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whe  re 


5Vp  is  the  incremental  plenum  volume  due  to  wave. 
Wp  is  the  width  of  the  plenum  chamber, 
dx  is  the  differential  length  of  the  plenum. 


The  total  incremental  volume  of  the  plenum  due  to  the  wave 
can  be  established  at  any  time  t  by  integration  over  the  plenum 


length,  L^.  Thus, 


5V  ( 

P 


/r2 

t)  =  W  A  /  sin(oj 
p  J  i 


(0  .  X  (0  .  X  ' 

t)  cos  +  cos  ( to  t)  sin  1 


V 


L 

— £ 
2 


w 


V 


w 


dx 


It  can  be  observed  that  the  first  term  on  the  right  hand  side 
of  the  equation  when  the  limits  are  applied  wiTl  yield  a  double¬ 
amplitude  sinusoid  (since  it  is  an  odd  function)  and  the  second 
term  will  yield  zero  (since  it  is  an  even  function) .  Thus, 


6Vp(t) 


sin(w  t) 
e 
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Now,  for  each  regular  wave  input  with  a  wave  length  of  L  and 

2ttV  IgL 

W  + 1  w 

an  incident  (circular)  frequency  of  =  — -  where  t^ie 


w 

coefficient  in  the  brackets  above  is  constant  and  thus  the  incre¬ 
mental  variation  of  plenum  volume  is  the  result  of  the  sinusoidal 
input,  5W  =  A  sin(a>et)  . 

The  transfer  function  of  the  incremental  plenum  volume,  Sv^ 
to  wave  height,  5w  becomes, 


sin 


where  the  coefficient  ki  (ok)  is  a  function  of  the  incident  frequency. 

Following  the  same  path  of  development,  the  transfer  function  of 
each  sidewall  volume,  5Vs  ,  to  wave  height,  SW  is. 


5V 

SW  =  k°(“i)  = 


'2W  V  \ 
s  w' 

V*i  ) 


sin 


where  W  =  width  of  each  sidewall, 
s 

The  transfer  function  of  the  total  incremental  sidewall  buoy¬ 
ancy  force,  6Fs  to  wave  height,  5W  is 


Then 


k0  (coi)  (2pwg) 


5F 

W 


sm 


k0  (<ai) 


2p  g 
w3 

M 


where  6F  =  6F  /M 
sm  s' 
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3.  System  Block  Diagram  and  Transfer  Function 

The  linear  heave  motion  equations  were  developed  in  Ref.  1  and 
are  repeated  here  for  convenience  in  Appendix  A.  The  block  diagram 
for  this  linear  system  is  shown  in  Figure  1  where  the  output  is  the 
incremental  variation  in  heave  (or  draft)  6z,  and  the  input  is  the 
sinusoidal  variation  in  wave  height  about  the  calm  water  level  of 
draft.  Note  that  a  positive  value  of  wave  height  reduces  the  effec¬ 
tive  plenum  volume,  5V^  ,  and  therefore  requires  a  negative  sign  in 
the  signal  transmission  as  shown  in  Figure  1. 

Also  it  should  be  noted  that  a  positive  value  of  wave  height 
increases  the  wetted  sidewall  volume,  5Vs  ,  which  cause  a  larger 
upward  force  and  therefore  requires  a  negative  sign  in  the  signal 

transmission.  The  effect  of  added  mass,  <SF  ,  was  also  included  in 

am 

this  mo^el.  However,  in  the  development  of  the  transfer  function  of 
C.G.  acceleration  to  wave  height  as  well  as  plenum  pressure  to  wave 
height,  it  was  shown  that  the  effect  of  added  mass  can  be  neglected, 
as  a  reasonable  first  approximation  and  therefore  added  mass  effects 
were  not  included  in  the  analysis  that  follows. 

The  algebraic  computation  of  the  wave  input  transfer  functions 
can  be  performed  either  by  block  diagram  manipulation  or  by  applica¬ 
tion  of  the  signal  flow  graph  method.  The  latter  approach  was  used 
with  the  details  contained  in  Appendix  B.  The  resulting  transfer 
functions  are  of  the  following  forms 

6 z  _  -  (a  s  +  b) _  (Draft/Wave  height)  (1) 

s 3  +  cs 2  +  ds  +  e 
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12. 


FIGURE  1.  BLOCK  DIAGRAM  OF  THE  LINEAR  HEAVE  EQUATION 

FOR  SINUSOIDAL  WAVE  INPUT 


/C.G.  acceleration^ 
\Wave  Height  )  {  ’ 


^  s3  +  cs2  +  ds  +  e 


<5z  _  -s2  (a  s  +  b) 


^Pb  _  s2 (a' s  +  b1 ) 


( 


Plenum  Pressure 

Wave  Height 


$W 


s 3  +  cs 2  +  ds  +  e 


The  coefficients  a,b,a'  and  b'  are  functions  of  the  incident 
wave  frequency  since  they  contain  the  k0  and  ki  terms  of  the  plenum 
volume-wave  height  transfer  function  and  sidewall  volume-wave  height 
transfer  function  developed  in  secion  III-A-2.  It  is  possible  to 
simplify  this  functional  relationship  in  k0  and  ki  by  a  convenient 
approximation  which  is  carried  out  below. 

From  Appendix  B, 


a  =  ko  +  kik2k3k4 


where  k 


k2  =  1/Vfa(0) 
k  3  =  YPb(0) 
k4  =  Ah/M 

For  the  XR-3  craft  operating  at  30  knots,  the  encounter  fre¬ 
quency,  cue  ,  can  be  approximated,  as  shown  in  Appendix  C,  over  the 
frequency  range  of  interest  as,  =  2<d^2.  Recognizing  that  = 
g/tiK  ,  it  is  possible  to  apply  these  substitutions  into  k0  and  ki 
to  yield 
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sin  aco 


k0 


(jW  \ _ 

\  M  /  am 


and 


ki  - 


■  <V 


sin  aco 
aco 


where 


A.  =  W  L 
b  p  p 


A  =  W  L 
P  s  P 


a  =  Lp/4g 
co  =  coe 


Observe  now  that  ko  and  ki  are  now  composed  of  constant  terms 
multiplied  by  the  sine  function  S^; A  plot  of  this  function 
over  the  encounter  frequency  range  of  interest  is  given  in  Figure 
2.  It  should  be  noted  that  the  magnitude  of  this  function  is  very 
nearly  unity  up  to  2  rps  then  attenuates  in  a  manner  that  approxi¬ 
mates  a  single  pole  at  8.4  radians/sec  up  to  15  rps  after  which  it 
drops  at  a  faster  rate  and  then  a  lobe  of  increased  magnitude 
occurs  with  a  peak  at  approximately  30  rps  and  repeated  attenuated 
lobes  at  multiples  of  the  plenum  length  as  reported  in  reference  7. 


The  coefficient 

’2AoP,.,9  YPb(0)Al 


a  = 


srw- 
M 


V,  (  0 ) M 
a 


and  in  like  manner 


sm  aco 
aco 


b  = 


2Aspwg  YPb(0)k6 


M 


#  Ijo T 

m 


sin  aco 
aco 


=  b'  = 


w°> 

V°) 


sin  aco 
aco 
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RANGE  OF  <j  VALUES  FOR  WHICH  THE  HEAVE  GAIN  APPROXIMATION  IS  VALID 
- MAXIMUM  RANGE - ►! 


g 


‘  > 


IS 


where  the  added  mass  coefficient,  kg  ,  that  originally  appeared  in 
coefficient  b  has  been  set  equal  to  zero. 

The  transfer  functions  now  become 


<5z  _  -a(s  +  b/a) 

3  2  . 

s  +  cs  •+  ds  +  e 


(lb) 


<Sz  _  -as  (s  +  b/a) 

6W  3  2  . 

s  +cs  +ds+e 


(2b) 


a1  s  (s  +  1) 

s3  +cs2  +ds+e 


(3b) 


where 


and 


a  = 


a  = 


b  = 


'2Aspwg  +  YPb(0)Ab 

M 

YPb<0)Ab 

M  J  aw 

’2Asp«g  Pb(0)k, 

M  VtQ) 

.  m 


sm  aco 


sin  aa) 
aco 


(since  yPb(0)Ab  >>  2AsPwg) 


sin  aco 
aco 


a1  = 


W0) 

vb(0) 


sin  aco 
aco 
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YPb(0) 

V  (0)  *• 
m 

2AsP„g  _  YP„«»V 
M  V.  ( 0 )  M 

D 


with  added-mass  coefficient 
k9  =  zero 


d  = 


^pb(0)V 

Vfc(0)M 


e  — 


YPb(0)k6 

v(0) 

m 


(since  yP,(0)A,  >> 


vb(0) 


fVw* 

.  M 


2V„g) 


The  value  of  b/a  is  dependent  upon  the  magnitude  of  the  air 
flow  rate  term,  k6  ,  thus 


b  .  2VwgYV0)ks  M 

a  Mvm<°>  ’  <2ASPW9  +  YPblO)^) 


which  will  reduce  by  the  above  approximation  to  yield, 


b 

a 


2Vwg 

W01 


Using  the  parameter  values  given  in  Section  IV,  the  transfer 
functions  are 


<5  z 

6W 


‘-3200.  (s  +  .045  k8) 
s3  +  11.3  kss2  +  2441.8s  +  145.3 


sin  ( .  155cj ) 
( .  155(o) 


6  z  _  2  <5  z 

<SW  ~  s  6W 


6Pb 
6  W 


2284  s2  (s  +  1) _ 

s3  +  11.3  k6  s2  +  2441.8s  +  145.3 


sin  (.155(o) 

(1.55co) 
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where 

a  = 

~2  (Ji)  2  (32.17)  +  1.4  (2140.8)  200  "1  sin  (4  -yy)  <u 

188.06  J  (5/32.17)  co 

b  = 

~2  (lr~)  2(32.17)  1.4(2140.8)  k6l 

4  0  sm  otco 

188.06(264.63)  aco 

a'  = 

200(1.  4)  (2140.8)  sin  aco 

262.44  J  aco 

c  = 

1.4(2140.8)  . 

264.63  J  6 

d  = 

2(—)  2(32.17)  H.420  (200)2 

188.06  188.06 

e  =  (11.326) (12.83)  k6 

Using  the  value  of  kg  =  6.28,  the  transfer  functions  become 


6  z 

6W 

sin  (  .  155co) 

-3200  (s  +  .28)  (  .  15 5co ) 

(s  +  .38)  (s  +  35.1  +  j  34.4)  (s  +  35.1  -  j  34.4) 

6  z 

6W 

2  6z 

3  Sw 

6Pb 

zw~ 

sin (. 155  ) 

2284  s2  (s  +1)  ( .155  ) 

(s  +  .38) (s  +  35.1  +  j  34.4) (s  +  35.1  -  j  34.4) 

where  the  factored  form  of  the  denominator  was  obtained  from  Section 
IV  A. 

Figure  3  shows  the  magnitude  spectrum  for  the  three  transfer 
functions  of  the  simplified  XR-3  model  for  the  value  of  k6  =  6.28. 
The  curve  for  acceleration  was  plotted  in  g  units  and  the  pressure 
curve  in  incremental  (gage)  pressure  values. 
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B.  Calculation  of  the  6  DOF  Model  Frequency  Response. 

Using  the  6  DOF  equations  for  the  XR-3  as  given  in  Ref.  5, 
the  output  response  to  a  small  amplitude  (10%  of  the  average  draft) 
sinusoidal  head-on  sea  was  obtained  on  the  IBM  360/67  computer. 

Two  types  of  wave  inputs  were  used;  a  single  frequency  input 
(regular  sea)  and  a  complex  frequency  input  (complex  sea) .  The 
band  of  encounter  frequencies  was  extended  beyond  that  used  in 
Ref.  6  to  include  the  range 

1  <  w  <  100  rps 
1.  Regular  Sea 

The  procedure  used  in  the  computation  of  output  response  to 
a  regular  sea  input  was  to  compute  35  seconds  with  the  adjustable 
thrust  option  to  arrive  at  a  steady  state  sinusoidal  response. 

The  average  value  of  thrust,  draft,  pitch  and  pressure  were  then 
used  as  input  for  the  operating  condition  of  the  final  run.  The 
steady-state  final  run  was  obtained  with  constant  thrust  at  a 
steady-state  speed  of  30  knots.  Data  was  taken  only  after  at 
least  30  seconds  of  compute  time  to  insure  true  steady-state 
conditions.  The  integration  method  used  was  a  f ixed-stepsize 
Runge  Kutta  which  yielded  about  50%  saving  in  compute  time  over 
the  adjustable  R-K  method  with  no  observable  loss  in  accuracy. 

The  ratio  of  draft  to  wave  height  amplitude  as  well  as 

C. G.  acceleration  and  plenum  pressure  to  wave  height  were  computed 
for  the  6  DOF  equations  with  the  results  plotted  in  Figure  4.  All 
3  frequency  response  curves  clearly  show  the  S1^  function 


20 


-+ - h 

O  1H9I3H  3AVM  o 

8  3anss38d 


8 


o 


H - h 

o  1H9GH  3AVM  0 

»  N0llVd31300V9  0  ~ 


H - 

w  1H9I3H  3AVM 

o  uvya 


2./ 


lobes  occuring  at  approximately  28  rps  and  above.  A  comparison 
of  Figure  4  and  Figure  3  clearly  shows  that  the  simplified  XR-3 
is  in  good  agreement  with  the  6  DOF  results.  For  frequencies 
below  28  rps,  the  transfer  functions  can  be  approximated  from 
the  data  of  Figure  4  using  the  straight  line  approximation  methods 
of  Bode  as  follows : 

For  draft,  the  frequency  response  follows  very  closely  the 
curve  of  Figure  2  for  the  sine  function  thus 


5  z  ^ 

m 


sin  aco 
aw 


where  a  =  Lp/4g 
w  =  we 


For  pressure,  the  response  is  approximated  by  the  transfer 
function 


!Ib 

5W 


sin 

ao i 


s 


For 

transfer 


C.G.  acceleration,  the  response  is  approximated  by  the 
function  (in  g  units)  as. 


Sz 

SW 


0.35 


sin  aw\ 
aw  J 


These  above  approximations  neglected  the  small  amplitude 
variations  that  occur  between  3  and  6  radians  per  second  that  are 
most  likely  introduced  by  the  pitch  mode  coupling  which  is  not 
accounted  for  in  the  simplified  heave-only  model  of  this  study. 


2 .  Complex  Sea 

The  complex  or  irregular  sea  conditions  were  approximated  by 
the  summation  of  10  regular  sea  components  as  provided  for  in  the 
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6  DOF  wave  subroutine  but  with  sufficient  amplitude  in  each  com¬ 
ponent  to  introduce  an  output  response  without  exceedinq  the  total 
height  limitations  for  maintaining  somewhat  linear  operation.  The 
transfer  functions  amplitude  ratios  were  computed  at  these  dis¬ 
crete  frequencies  using  the  Fast  Fourier  Transform  of  Reference  4 
and  the  values  plotted  in  Figure  4  as  square  points.  It  can  be 
seen  that  each  of  the  discrete  frequencies  produces  values  that 
are  in  reasonable  agreement  with  the  regular  sea  input  results.  This 
also  has  been  confirmed  by  Booth  in  Reference  4.  It  should  be 
noted  that  the  complex  sea  input  frequency  response  computation  is 
a  single  run  technique  but  limited  by  the  following 

1)  Discrete  frequencies  must  by  necessity  be  limited 
in  number  and  so  chosen  that  they  are  harmonics  of 
the  fundamental  frequency. 

2)  Each  frequency  must  introduce  a  reasonable  ampli¬ 
tude  to  produce  a  response  which  in  total  must  be 
limited  to  maximum  composite  wave  height. 

3)  Steady-state  speed  must  exist  since  encounter  fre¬ 
quency  must  be  a  multiple  of  the  fundamental  fre¬ 
quency  to  reduce  the  effect  of  leakage  (a  problem 
with  the  Discrete  Fourier  Transform  method) . 
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IV.  PLENUM  AIR  FLOW  RATE  PARAMETER  STUDY 


During  the  study  conducted  on  scaling  the  heave  motion 
equation  as  reported  in  Reference  2,  it  was  observed  that 
the  air  flow  influx  as  well  as  the  leakage  rate  determined 
the  characteristic  response  of  heave  motion  for  a  step  weight 
disturbance.  It  is  of  interest,  therefore,  to  study  the 
sensitivity  of  the  craft  dynamics  to  these  parameters  (air  flow 
influx  and  leakage  rate) .  The  characteristic  equation  of 
our  simplified  linear  model  is  a  function  of  the  parameters,  so 
a  root  locus  study  permits  analysis  of  their  effect  on  root 
locations.  Alternately  the  transfer  functions  may  be  recomputed 
and  the  frequency  responses  plotted  and  analyzed.  Finally  the 
six  degree  of  freedom  simulation  program  may  be  exercised  to 
provide  both  frequency  response  data  and  step  weight  removal 
transients . 

A.  Characteristic  Equation  Root  Locus 

The  linearized  equations  (see  Appendix  A)  yield  the  following 
characteristic  equation 

s3  +cs2  +ds+e=0 

which  in  terms  of  the  system  parameters  becomes, 
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represent  the  adjustable 


parameter  for  the  characteristic  equation  which  can  be  re-organized 
into  the  root  locus  form  for  a  ir-type  locus , 


K, 


(1.4) (2140.8) 
264.63 


_  2  ,  2(1.99)  (32.17)  (75)  (32.17) 
5  (4)  (6050) 


-1  = 


\ 


s<s‘  + 


2  (75)  (32.17)  (1.99)  (32.17)  ,  (  1. 4 )  ( 200 )  2  ( 2140 . 8)  (  32 

4(6050)  (262.44)6050 


nl 


Then  the  equation  for  the  root  locus  becomes 


K. (11.33) (s2  +  12.83) 

— -  =  -1 

s (s2  +  2441. 8) 


To  carry  out  the  root  locus  study  numerical  values  are 
required  for  the  parameter.  We  obtained  the  values  used  in 
the  following  study  from  the  six  degree-of-f reedom  program. 
These  values  were  also  used  in  the  simulation  studies  for 
the  10%  weight  removed  transient;  thus  the  root  locations 
can  be  correlated  with  the  transient  response. 
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An  IBM  360/67  subroutine  called  ROOTLO  was  used  to  compute 
and  plot  the  locus  of  roots  with  the  results  shown  in  Fig.  5  and 
in  Table  I.  The  values  of  used  in  this  study  covered  the 
range 

6  <  K£  <  22 

with  the  root  locus  of  Fig.  5  covering  the  range 

5  <  Kz  <  50 

From  Fig.  5,  and  Table  I,  it  can  be  seen  that  the  roots  of 
the  characteristic  equation  for  small  values  of  K,  are  one  real 
pole  near  the  origin  and  a  complex  pair  far  removed  from  the  origin. 
For  larger  values  of  K,,  the  complex  pair  become  real  roots  with 
one  root  approaching  the  root  near  the  origin  and  the  other  moving 
towards  negative  infinity.  Note  in  particular  that  the  real  root 
near  the  origin  moves  a  relatively  short  distance  during  the  one 
decade  change  of  the  parameter  K^. 

In  summary,  it  is  noted  that  the  real  root  near  the  origin 
is  not  very  sensitive  to  K  changes  in  this  range  of  values , 
whereas  the  two  other  roots  are  quite  sensitive  to  changes  in  K^. 

Since  K „  is  a  function  of  air  flow  rate  influx  (nk,)  and  air  flow 

CnA*  FI 

rate  leakage  -  .1- - — r-  ,  it  is  clear  that  the  oscillatory 

p  a  V2Pb(0) 

modes  in  the  ship  dynamics  depend  on  these  flow  rates. 


26 


CO 

o 

o 


X 


+ 


+ 


o 

o 


+ 


+ 


+ 


+ 


CM 

o 

o 


+ 


1 -008  ^ 


-006 


-ool;  1  1 


L^Hn|||||A|00 


+ 


O 

o 


+ 


+ 


+ 


+ 


+ 


o 

o 


002 


004 


006 


X 


CO 

o 

o 


GO 

o 

o 


X-SCALE=2 . 00E+0 1  UNITS  INCH. 

Y-SCALE=2 . 00E  +  0  l  UNITS  INCH. 

BOX  70  GERBfl  PARAMETER  STURY 
FIG. 5  ROOT  LOCUS  FOR  KL  ADJUSTABLE 

27 


Table  I  -  Root  Location  for  Various  Values  for  the  Characteristic 
Equation,  (s  +  a)  (s  +  6)  (s  +  y)  =  0 


a 

6 

y 

5.00 

0 

+ 

j 

49.4 

0 

_ 

j 

49.4 

0.0 

6.28 

-35.1 

+ 

j 

34.4 

-35.1 

_ 

j 

34.4 

-0 .38 

7.88 

-44.0 

+ 

j 

21.5 

-44.0 

_ 

j 

21.5 

-0.47 

9.88 

-81.3 

+ 

j 

0.0 

-29.2 

j 

0.0 

-0.60 

12.4 

-119. 

+ 

j 

0.0 

-19.6 

j 

0.0 

-0.77 

15.6 

-160. 

4* 

j 

0.0 

-14.2 

j 

0.0 

-0.99 

19.5 

-208. 

+ 

j 

0.0 

-10.4 

i 

0.0 

-1.31 

24.5 

-267 

+ 

j 

0.0 

-7.3 

— , 

j 

0.0 

-1.  82 

28.5 

-313 

+ 

j 

0.0 

-5.3 

j 

0.0 

-2.50 

33.2 

-366 

+ 

j 

0.0 

-3.31 

+ 

j 

1.  44 

-3.31 

-  j  1.44 

38.6 

-428 

+ 

j 

0.0 

-2.84 

+ 

j 

2.23 

-2.84 

-  i  2.23 

44.9 

-500 

+ 

0.0 

-2.43 

+ 

2.66 

-2.43 

-  3  2.66 
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B.  Sinusoidal  Sea  Input  Frequency  Response 

Both  the  simplified  linear  model  and  the  six  degree  of  freedom 
program  were  used  in  this  part  of  the  study.  Figure  6  shows  the 
frequency  response  of  the  linear  model  for  high  leakage  conditions. 
Comparison  of  Figure  6  with  Figure  3  shows  one  effect  of  the 
leakage  conditions:  with  high  leakage  the  higher  frequency 
motions  are  substantially  attenuated.  Note  that  for  encounter 
frequencies  less  than  8  rad/sec  (wg  <  8)  the  craft's  vertical 
motion  is  not  at  all  sensitive  to  air  flow  and  leakage  rates, 
but  for  wg  >  10  the  CG  acceleration,  pressure  and  heave 
displacement  are  all  substantially  reduced  for  the  case  of  increased 
leakage . 

Figure  7  shows  the  frequency  response  of  the  six  degree  of 
freedom  model  for  high  leakage.  Comparison  of  Figure  7  with  Figure  4 
shows  the  same  general  characteristics  observed  with  the  linear 
model:  the  response  to  frequencies  above  10  radians  per  second 

are  significantly  attenuated  when  the  air  flow  is  increased. 
Frequencies  below  10  radians  per  second  produce  somewhat  larger 
responses  when  the  airflow  is  increased.  We  believe  this  is  a 
result  of  pitch  motions,  as  discussed  in  the  following  paragraphs. 

C.  Comparison  of  the  Linear  Model  with  the  6DOF  Model  Using  the 

Frequency  Response  Plots 

One  purpose  of  this  study  is  to  assess  the  value  of  the  linear 
model  as  a  tool  for  studying  the  CAB  and  predicting  the  salient 
features  in  its  behavior.  Alternately,  comparison  of  results 
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between  the  linear  model  and  the  6  DOF  model  (or  experimental  data 
when  available)  should  lead  to  an  understanding  of  phenomena 
which  are  peculiar  to  the  CAB.  On  Figure  8  the  d.g.  acceleration 
curve  of  the  simplified  model  and  that  of  the  six  degree  of 
freedom  model  are  plotted,  using  the  data  for  K,  =  19.99. 

At  low  frequencies  the  six  degree  of  freedom  curve  indicates 
much  larger  accelerations  than  the  linear  model.  At  high 
frequencies  both  models  are  in  good  agreement  (the  null  frequencies 
are  not  exactly  the  same,  indicating  that  the  effective  plenum 
length  is  not  precisely  the  geometric  length) .  In  the  low 
frequency  range  we  note  that  the  six  degree  of  freedom  frequency 
response  has  a  slight  "dip"  or  anti -resonance  for  4  <  w  <  6, 
while  the  linear  model  does  not.  This  frequency  range  includes  the 
natural  frequency  of  the  pitch  mode,  which  is  included  in  the  six 
degree  of  freedom  model,  but  not  in  the  linear  model.  We  know 
that  the  sinusoidal  sea  will  excite  the  pitch  mode,  so  we  are 
led  to  a  plausible  explanation  for  the  discrepancy  between  the 
results  at  low  frequencies.  It  is  simply  this  -  with  sinusoidal 
excitation  in  the  range  of  frequencies  about  pitch  resonance, 
the  pitching  motion  will  have  large  enough  amplitude  to  increase 
the  air  leakage.  With  the  bubble  volume  thus  reduced  the 
incompressible  wave  will  cause  higher  air  pressures  and  higher 
heave  acceleration.  At  high  frequencies,  above  the  natural 
frequency  in  pitch,  leakage  is  not  increased  and  the  linear  model 
is  a  reasonably  good  one. 
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We  also  note  that  at  low  frequencies  the  slopes  of  the 
curves  on  Figure  8  are  not  the  same.  This  can  be  attributed 
to  variation  in  the  leakage  as  a  function  of  frequency. 

D.  Simulation  of  Transient  Response  due  to  Step  Removal  of  Weight 
The  removal  of  weight  by  a  step  change  is  equivalent  to  a 
negative  step  change  in  acceleration.  The  transfer  function  for 
draft  pressure  and  C.G.  acceleration  for  the  step  weight  removal 
can  readily  be  found  by  Mason's  Gain  Rule  using  the  signal  flow 
graph  given  in  Figure  B-l  of  Appendix  B.  The  transfer  function 
for  draft  becomes 

Apb<°)V 
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5a  D(s) 


where  6z  =  incremental  draft 
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These  transfer  functions  are  for  the  linear  simplified 
model,  for  which  we  have  the  root  locations  (factors  of  D(s))  in 
Table  I  with  corresponding  root  locus  on  Figure  5.  The  transfer 
functions  also  indicate  the  presence  of  zeros ,  some  of  which 
are  functions  of  the  leakage,  K  . 

In  studying  the  weight  removal  behavior,  simulation  was  done 
with  the  six  degree  of  freedom  model,  and  the  results  are 
conveniently  analyzed  if  we  consider  Table  I  and  Figure  5 , 
remembering  also  that  the  simulation  model  contains  additional 
pitch  mode  effects.  Figure  9  is  a  set  of  16  curves  covering  the 
range  of  air  flow  conditions  from  6.1  <  <21.99,  but  recorded 

for  less  than  a  1.0  second  duration  after  removal  of  the  weight. 

This  expands  the  scale  of  the  initial  transient,  and  permits 
correlation  with  Table  I  and  Figure  5.  Figure  10  is  a  set  of 
12  curves  covering  the  same  air  flow  conditions,  but  recorded  for 
a  5.0  second  duration,  which  permits  us  to  recognize  the  effect 
of  pitch  coupling. 

Consider  first  the  c.g.  acceleration  curves  of  Figure  9 
(Al,  A2,  A3,  A4 ) ,  noting  from  Table  I  that  for  each  of  these 
figures  the  value  results  in  a  small  negative  real  root,  but 
gives  complex  roots  only  for  Figure  9-A1,  where  c  =  0.7  and  wn  >  50. 
Each  of  these  figures  shows  a  very  fast  initial  transient,  there 
is  no  apparent  contribution  from  the  small  real  root  (Table  I, 
column  y) ,  which  may  have  been  cancelled  by  the  zero  of  the 
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transfer  function.  Comparison  of  these  four  initial  transients 
shows  very  good  agreement  with  the  root  values  of  Table  I  for 
low  leakage  (Figure  9-A1) .  The  initial  transient  shows  a  small 
overshoot,  about  as  expected  for  roots  with  ?  -  0.7.  As  K, 
increases  (Figures  9-A2,  A3,  A4)  the  transient  does  not  overshoot 
(roots  are  all  real)  and  the  settling  time  increases  which 
agrees  with  the  change  in  root  values  in  column  8  of  Table  I. 

Continuing  the  analysis  of  Figure  9,  the  z-displacement  and 
plenum  pressure  curves  show  an  approximately  exponential 
variation  which  is  very  slow  but  speeds  up  with  increased  leakage. 
This  agrees  qualitatively  with  the  increase  in  value  of  the 
smallest  real  root  (column  y,  Table  I).  The  pitch  angle  variation 
of  Figure  9  shows  a  substantial  variation,  but  the  short 
duration  does  not  supply  enough  information  for  further  analysis. 
We  therefore  proceed  to  Figure  10. 

Let  us  start  the  discussion  of  Figure  10  with  the  pitch 
angle  response  curves,  10-Dl,  D2 ,  D3.  On  all  of  these  curves  it 
is  clear  that  the  pitch  transient  is  an  exponentially  damped 
sinusoid  plus  another  exponential  response.  The  frequency  of 
this  pitch  mode  is  approximately  5  rad/sec,  and  increases  slightly 
with  increased  K^.  We  cannot  correlate  this  with  the  linear 
model,  for  it  does  not  include  the  pitch  mode.  The  second 
exponential  component,  however,  appears  to  be  due  to  the  small 
real  root  of  the  heave  mode  (Table  I,  column  y)  and  its  time 
constant  seems  to  decrease  with  increasing  ,  which  agrees  with 
the  root  variation  in  Table  I. 
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From  Figure  10  it  is  also  clear  that  the  pitch  mode  couples 
into  c.g.  acceleration,  z-displacement ,  and  plenum  pressure. 

This  coupling  increases  with  K^.  The  small  real  root  is  also 
seen  to  contribute  —  it  is  obviously  dominant  in  the  z-displacement 
and  pressure  curves,  and  can  be  recognized  in  the  c.g.  acceleration 
curve  of  Figure  10-A3. 
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V.  CONCLUSIONS 


A  few  specific  conclusions  may  be  drawn  from  the 
results  of  this  study: 

1.  The  major  effect  of  increasing  air  flow  rates  is 

to  attenuate  high  frequency  components  in  the  response 
of  the  craft.  This  is  clearly  shown  by  the  high 
frequency  response  of  both  the  linear  and  6  D.O.F. 
models.  It  is  also  shown  by  the  step  weight  removal 
transients.  Therefore,  to  reduce  the  vibrations  due 
to  vertical  accelerations  at  high  encounter  frequencies, 
the  air  flow  rate  should  be  maintained  as  high  as 
practical. 

2.  The  pitch  mode  couples  strongly  into  the  vertical 
motions  over  the  encounter  frequency  range  of  1  to  10 
rad/sec.  as  shown  by  the  6  D.O.F.  model.  Therefore, 
pitch  motion  should  be  added  to  the  linear  model  to 
account  for  the  changes  in  air  flow  leakage  rate 

due  to  the  variations  in  the  pitch  angle. 

3.  The  linear  model  (Heave  motion  only)  is  sufficiently 
accurate  for  many  studies  and  should  be  used  where 
applicable  to  aid  in  the  analysis  and  design  of  CAB 
craft  characteristic. 

4.  The  linear  model  reduces  both  the  CPU  execution  time 
and  the  operating  class  in  the  job  queue  resulting  in  a 
reduce  interaction  (turnaround)  time. 
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5.  The  complex  sea  method  of  developing  the  frequency 
response  of  the  6  DOF  model  using  the  Fast  Fourier 
Transform  provides  approximate  results  that  can  be 
quite  inaccurate  unless  special  care  is  exercised 
in  its  application. 
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VI.  RECOMMENDATIONS 


One  of  the  long  range  goals  of  the  analytical  studies 
of  the  CAB  craft  at  NPS  is  to  develop  a  simplified  model 
that  will  provide  a  reasonable  computer  execution  time  for 
control  system  design  and  evaluation  in  the  complex  sea 
operating  environment.  We  feel  that  we  have  a  reasonable 
start  in  this  direction  and  have  the  following  recommendation 
for  future  studies : 

1.  The  pitch  mode  should  be  modelled  and  incorporated 
in  the  Linear  Heave-Only  Model.  Mr.  McIntyre  in 
Reference  7  has  made  the  first  effort  in  this 
direction  and  a  current  study  by  Mr.  Barnes  (Reference  8) 
has  improved  on  the  model  and  shows  promise  of 
providing  the  additional  information  missing  from 

the  Heave-Only  Model. 

2.  Upon  the  completion  of  the  Linear  Pitch-Heave  Model 
and  its  regular  sea  frequency  response  validation 
with  the  6  D.O.F.  model,  the  complex  sea  frequency 
response  should  be  computed  and  compared  for  both 
accuracy  and  speed  of  computation. 

3.  Continue  experimentation  with  both  software  and  hardware 
techniques  for  the  analysis  of  the  complex  sea  results 
both  for  the  6  D.O.F.  model  and  the  actual  XR-3  craft 

to  provide  improved  results  of  the  complex  sea  frequency 
response  analysis. 
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4. 


Develop  deterministic  input  tests  for  the  XR-3 
craft  for  validation  of  the  SESCOMP  model.  The 
weight  removal  transient  is  one  such  test  that 
shows  promise  for  validation  of  the  Pitch-Heave 
mode  equations. 

5.  Investigate  effect  of  the  addition  of  a  membrane 
in  the  plenum  chamber  to  reduce  the  variations  in 
the  vertical  motions  of  the  XR-3  craft.  The 
actual  XR-3  craft  as  well  as  the  6  D.O.F.  program 
and  the  linear  pitch-heave  model  would  be  used 
in  this  study. 
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APPENDIX  A.  DESCRIPTION  OF  THE  SIMPLIFIED  XR-3  MODEL 


A.  Nonlinear  System 

In  this  model  of  the  XR-3,  pitch  variations  are  reduced 
to  zero  by  assuming  1)  the  center  of  pressure  (C.P.)  directly 
under  the  center  of  gravity  (C.G.),  2)  the  sidewall  symmetrical 
about  the  C.G.  with  uniform  rectangular  cross-section  from 
bow  to  stern  and  3)  the  pitch  moments  of  the  seals  and 
aerodynamics  cancel  each  other.  The  lift  force  of  plenum 
pressure  and  sidewall  buoyancy  oppose  the  craft  weight.  All 
other  lift  and  drag  forces  are  neglected  since  constant  speed 
conditions  are  assumed. 

In  addition  it  will  be  assumed  that  the  rear  seal  maintains 
a  constant  leakage  area,  that  is,  the  seal  follows  the  water 
level  at  a  fixed  separation. 

Figure  A1  shows  a  top  and  sideview  of  the  simulated  CAB 
craft.  All  dimensions  are  chosen  to  approximate  those  of 
the  XR-3  craft  in  order  to  compare  the  results  of  this 
analysis  to  the  L  &  M  program  solution  under  similar  operating 
conditions. 

The  equations  of  motion  for  this  system  are  given  below. 


Orifice  Leakage  Rate,  <3out  -  cnAz 


cu. ft. 
sec . 


(1) 


Fan  Map  Input  Rate ,  q . 


in 


n 


cu. ft . 


sec . 


(2) 
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FIGURE  Al.-Top  and  Side  View  of  Simplified  XR 


-  -  -  " 

— 

—  -  . 

Absolute  Plenum  Pressure,  P 

b  =  pa  Psf 

(3) 

(Adiabatic  Process) 

b  lVa) 

Plenum  Volume,  V,  =  ( V  -  A, 

d  no 

l cu.ft. 

(4) 

Plenum  Air  Flow  Rate, 

M.  =  p 
b  a 

(a  a  lbs’ 

^qin  qout)  sec. 

(5) 

Heave  Acceleration,  Z 

\M 

!e  -  Ft\  ft. 

M  M  /  sec . 2 

6) 

Heave  Velocity,  Z  =  j 

[fi  at 

1  ft. 

1  sec . 

(7) 

Heave ,  Z  = 

\fi  at 

+  Z (o)  ft. 

(8) 

L“o 

J 

P lenurrt  Pressure  Lift 

Force , 

F  =  A,  P,  lbs. 

p  b  b 

(9) 

Buoyancy  Lift  Force, 

II 

Cm 

2(A4-fd)pgJ  lbs. 

(10) 

Plenum  Gage  Pressure, 

,  P.  =  (PK  -  P  )  Psf. 

b  b  a 

(ID 

Draft, 

*d  = 

Z  +  Zs  feet 

(12) 

The  system  parameters  and  constants  are  listed  below. 


Adiabatic  process  coefficient. 
Leakage  area. 

Leakage  orifice  coefficient, 
Air  density. 

Atmospheric  pressure. 


Y  =  1.4 


A^  =  Yl  ZL  =  0.438 

C  =0.90 
n 

p  =  .002378 

cl 

P  =  2116. 

cl 
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Plenum  area. 

Empty  plenum  volume. 

Craft  weight. 

Keel  line  area  of  sidewall, 
Draft,  initial  condition. 
Density  of  water, 

Gravity  acceleration. 

Center  gravity  location. 
Number  of  fans. 

Steady  state  fan  output. 
Fan  map  slope. 


Afa  =  200. 

V  =  383. 
n 

W  =  6007 

A  =  75/4 

i.(o)  =  0.5 
d 

p  =  2 . 

g  =  32. 

Zs  =  2.5 
n  =  5. 


Qi(o)  =  35. 
kq  =  1.0 


The  schematic  diagram  of  signal  flow  is  shown  in  Figure  A2 . 
The  coordinate  system  and  initial  conditions  for  heave  and 
draft  are  shown  in  Figure  A3 . 


/|\ 


Zs  =  2.5  f t . 


7\ 
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_  . 
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£d(o)  =  0.5ft. 
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water  level 
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Figure  A3  -  Coordinate  System  and  Initial  Condition 
of  Heave  and  Draft. 
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FIGURE  A2 . —  Schematic  Diagram  of  Signal  Flow. 


It  should  be  noted  that  the  range  of  heave  motion  is 

-2.5  <  Z  <  -0.5  since  -2.5  feet  would  put  the  water  level  at 

the  keel  line  and  -0.5  feet  would  put  the  water  level  at  the 

plenum  roof.  Draft  is  measured  from  the  keel  and  would  have 

the  range  0  <  <  2.0.  Also  note  that  if  the  craft  "drops" 

into  the  water,  increases  and  Z  becomes  less  negative,  that 

is,  Z  increases  in  the  positive  (downward)  direction. 

At  steady-state  condition,  the  plenum  pressure  assumed 

is  P,(o)  -  24.  PSF  producing  a  F  (o)  =  4800  lbs.  For  the 
d  p 

draft  of  £3(0)  =  0.5  feet,  the  buoyancy  force  F  (o)  =  1200  lbs. 

The  air  mass  in  the  plenum  for  the  above  steady-state 

where  q.  (o)  =  q  ,  (o)  =  55  cu.  ft. /sec.  is  M,  (o)  =  0.71878 
in  out  o 

pounds . 

B.  Linear  System 

The  linear  system  equations  are  arrived  at  by 
application  of  the  Taylor  Series  expansion ^  about  the 
steady-state  operating  point.  The  development  is  shown  below 
for  the  two  differential  equations. 


a)  Steady  State  and  incremental  variations  are  noted  as  follows 
using  heave  as  an  example. 

z  =  z  (o)  +  6 z 
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Equation  (1) ,  let  W/M  =  Wm;  Fp/M  =  Fpm 

and  F£/M  =  F£m 


then, 


z<°>  +  -  »/«  -  <Fpm(0)  +  5Fpm)  -  (Fta(0)  +  5F«m> 


or. 


sz  *  -5Fpm  '  5Flm 


(la) 


Equation  (2) ,  let  m/p^  =  vm 


then , 

V  (o)  +  <5V  =  nQ.  -  nk  (P,(0)  +  <SP.) 

m  *  n\  -  4  D  D 


C  A  2Pb(0)  55 
nt  (—^ - ) 


2P.(o)  _k 
C  A  r,  H  (— - )  * 

n  c  p 


(— )  5PV 


or. 


•  _  C  An  ^  a.  k  _ 

6V  =  -  (nk  )  <SP. - -  ( -  )  ^  6P 

m  q  b  °a  2P,  (o)  b 


(2a) 


Next,  the  auxiliary  algebraic  equations  are  linearized. 

V°>  -  5FP  -  Ab  (Fb  <0)  -  SFb' 


«F„  ■  A.  «Pk 

p  b  b 


(3a) 
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F£  (o)  +  6FC  =  2A6  Pwg  (*d(o)  +  &ld) 

6F£  =  {2A5  pwg)  Ud  (4a)  | 


Eq. (5a)  and  (6a)  have  already  been  developed  they  are: 


6q.  =  -(nk  )  6P, 

^in  q  b 

(5a)  | 

o  .  _ 

6q  .  =  -  -  ( - - -  p  6P 

out  0  -)n  b 

a  2P  (°) 

! 

(6a  ! 

i 

1 

For  equation  7a 

(Pb(o)  +  5Pb) 

P 

a 


<5P,  =  P  y 

b  a 


V°> 

Vb(o) 


6  V 


m 


Vm(o) 


V°> 

vb(o) 


6  V, 


Vb(o) 


and  equation  8,  9,  and  10  follow. 
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Vb(o)  +  6Vb  =  Vn  -  Ab  (£d(o)  +  Ud) 


6vb 

(Afe)  6-^d 

(8a) 

(o)  + 

5Fb  =  pb(o>  +  «pb  -  pa 

6pb 

=  6Pb 

(9a) 

(o)  + 

Si,  =  z(o)  +  St  +  Z, 
a  4 

Si, 

d 

«  5  z 

(10a) 

l 

Combining  (la),  (3a),  (4a),  (7a)  and  (8a)  as  follows: 


-  -  ft) W 


S  z 


recall  that  Si,  =  S z 
d 


St.  -  -  b  \  y  P.  (o) 

—  i  b 


[vm(°) 

L  m 


6  V  - 
m 


-L_  (-A.  )  «zT 


vb(o) 


'  (2  t)  ** 


M 


6  z  =  - 


M 


2Vw^  +  (rPb(o)Ab2V 

/  \  MV,  (o)  /. 


MV,  (o) 

D 


6  z  + 


~'Pb(°)Abl  6V 

MVm(o)  J 


m 


Also  for  equation  (2a) 


6  V  = 
m 


nk^  +  ^n^t  /  pa  \ 

% 

ypk(o) 

6  V 

m 

-  6Vb 

Pa  \2P,  (o)/ 

1 

D 

_V  (o) 
m 

v^srJ 

6  V  = 
m 


nkq  +  Cn Rl  ( ea  ¥  yPh(o)  (  1  \ 

~  \^rzJ  j  D  L\  vm(o)/ 


a  2Pb (o) 


6  V 


m 


+  /  Ab  \  6z 

\V°>/ 


and  for  equation  (7a) 


5Pb  =  YPb(o) 


1  6  V  +  b  6  z 

m.  vTToT 

m  d 


Now  define  the  state  vector 

T 


x  A  [  6  z  6  z  6  V 


m 


and  output  variable,  y  =  6Pb  with 


y  A  Cn  X!  +  c13  x  3 


Then  the  state  equations  are 
xi  =  x2 

x2  =  a21  x3  +  a2  3  x3 
x3  =  a 3  j  x3  +  a 3  3  x3 
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where 


c,  =  ’fPb<°)Ab 
Vb(o) 


and 


_  _  yPK(o) 

C 1  3  -  b.  . 

V_  (o) 


m 


a2  1  “ 


/2A,p 

-[  AW 
V  M 


g  +  YPb(o)Ab: 


MVb(0) 


■) 


a2  3  ~ 


=  —  Apb(0,M 
\  “'.I"'  / 


a3  1  “  “ 


YPb (o) 


nk  +  CnA£ 

q  -7 — 


a  \J 


2Pb(o) 


*b 


a  3  3  ■  -  Ypb (°) 


i  .  C  A  p 
nk  +  n  £ 

q 


a  \| 


2Pb(0) 


Vo) 


1 

V  (o) 

m 


The  characteristic  equation  for  this  system  is  given  by 


(4 1  -  A)  where 


o 

1 

o 

A  = 

a  2  1 

o 

a  2  3 

a  3  1 

o 

a  3  3 

Then  (A I  • 

-  A)  = 

A 

-1 

o 

~a2  1 

<s 

“a2  3 

_“a  3  1 

o 

4  -  a33- 

and  the  characteristic  equation  is 

3  2 

A  -  a33  A  -  a2i  4  +  (^2i  a33  -  a23  a31)  =  o 

Using  the  given  numerical  values  the  characteristic  equation  is 

3  2 

A  +  62. A  +  2183A  +  789  =  o 

This  cubic  equation  will  factor  into  one  real  root  near  the 
origin  and  a  complex  pair  far  removed  from  the  jw  axis.  (See  Sec. 
IVA)  Thus  a  step  weight  response  will  have  a  complex  pair  of  poles 
which  can  be  approximated  by  the  quadratic  form, 
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4  —  3-3  3  ^  ~  ci  2  1  “  O 

For  the  given  numerical  values  it  can  be  shown  that  a 
rough  approximation  for  the  coefficients  are: 


-a 3  3  =  nk  y  P,(o)/V  (o) 
”  eq  b  m 


-a21  =  _b_  Y  Pb(o)/V  (o) 
M 


Using  the  above  approximation  yields  a  natural  frequency 


of : 


Wn  =  ~  = 


Ab  Y  Pb(o)/V  (o) 
M 


w  =  K  , 
n  n  \ 


Pb(o)/Vb(o) 


(13) 


where  Kn  =  Ab  \J  y/M 


and  a  damping  factor 


j.  =  _a3  3 

2  w 
^n 


n"kq"  y2  Pb2 (o)  /  V2m (o) 

4  Y  Pb(o)  /  Vb(o) 

M 


e  =  k4 


Pb(o)  Vb(o) 


V*  (o) 

m 


(14) 


where 


I - 

f  ■  2aT  \i  Ym 
b 


K 
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APPENDIX  B.  COMPUTATION  OF  THE  WAVE  INPUT  TRANSFER  FUNCTION 


The  block  diagram  of  Figure  1  given  in  Section  IIIA  can  be 
drawn  as  a  signal  flow  graph  as  shown  in  Figure  B-l.  To  obtain 
the  system  transfer  functions,  the  Mason  Gain  Rule  (MGR)  can  be 
applied  as  follows: 

A  —  1  —  ( L  i  +  1>2  +  Lg)  +  L2L3  +  Li  2  L  4 

where 

Li  =  -k2  k3  k4  k8  A-2 

L2  =  ~k3  k6  k7  a"1 

L3  =  -k5  A-2 

Lu  =  -kq  A  ~ 1 


For 

the  incremental  heave  (draft) , 

6z,  the  MGR  yields 

for  the 

wave  input ,  6W 

6z  =  i-  (Pi  Ai  +  P2A2)  6W 

where 

P\  =  -k0A~  2 

A 1  —  1  ■“  L2 

thus , 

P2  =  -ki  k2  k3  k4  A**2 

A  2  =  1 

6z 

-  I  [(-ko^'2>  (1  + 

+  (-k 1k2k3k4A~2 )  (1)] 

6W 


where  the  gain  coefficients  for  any  particular  operating 
conditions  are: 


k.  sin  f'J&W 


2  V 


w 


M 
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9H 


1 


ki 

= 

k2 

= 

^b 

k3 

= 

^Pb 

k4 

3 

VM 

k5 

= 

2V, 

/ 

k6 

= 

(^q 

k7 

3 

k8 

= 

> 

k9 

f  UCvp, 

\  8M 

w 


b  (o) 


sin 


2V, 


w 


pa 


b  (o) 


By  applying  substitution  and  reduction,  the  transfer  function 
becomes , 


6z  _  -(aa  +  b) _ 

<$w  J"5  +  c4 2  +  da  +  e 


with 

a  =  (k0  +  ki  k2  k3  k4) 

b  =  k0  k3  k6  k7 

c  =  (k  3  k6  k7  +  k9) 

d  =  (ks  +  k2  k3  k4  kg  +  k3  k6  k7  k9) 

e  =  k 3  k s  k6  k7 
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It  should  be  noted  that  the  transfer  function  of  draft 
to  wave  height  has  a  numerator  whose  coefficients  "a"  and  "b" 
are  functions  of  the  incident  wave,  and  therefore  will 

change  in  value,  for  each  frequency  of  the  sinusoidal  sea  input. 

In  addition,  it  should  be  noted  that  the  incremental  heave 
(draft)  is  180  degrees  out  of  phase  with  the  wave  input  as 
would  be  expected  since  a  positive  value  of  wave  height  produces 
(1)  increased  buoyancy  and  (2)  reduction  in  plenum  volume  which 
increases  the  plenum  pressure  and  thus  the  lift  force. 

Further,  it  is  observed  that  the  denominator  polynomial 
is  of  third  order  as  would  be  expected  since  the  characteristic 
roots  of  this  system  are  obtained  by  solution  of  the  equation 

•6  3  +  c6 2  +  d4  +  e  =  0 

as  given  in  Section  IV-A.  The  transfer  function  of  c.g. 
acceleration,  <5z,  to  wave  height,  6w,  is  obtained  from  the 
signal  flow  graph  as  follows. 


6z 

<5W 


-k1k2k3k4(l)  -  k0(l  +  k3k6k74-1) 


A 


•;  -t>  2  ( a-4  +  b ) 

0  z  _  _ 

-5  3  +  c4 2  +  d4  +  e 

where  the  coefficient  a,  b,  c,  d  and  e  are  the  same  as  those  in 
the  draft-wave  height  transfer  function. 
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The  minus  sign  appearing  in  this  transfer  function  represents 
the  180  phase  shift  which  is  also  shown  in  the  draft-wave  height 
transfer  function  as  would  be  expected. 

In  like  manner,  the  transfer  function  of  plenum  pressure, 

6Pb ,  to  wave  height,  6W  is 

6Pb  kik2k3(l  +  k54-2  +  k94_1)  -  k0k2k3k84”2 


6Pb  _  (a 1 4 2  +  bU  +  c1)  =  i2(aU  +  b1) 

A  A 

ki  k2  k3 
ki  k2  k3  k9 

k2  k3  (kiks  -  kokg)  =  zero 


where 


6W 
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APPENDIX  C.  RELATIONSHIP  OF  ENCOUNTER 
FREQUENCY  TO  INCIDENT  WAVE  FREQUENCY 
FOR  THE  XR-3  CRAFT 

Booth  (in  Appendix  A  of  Reference  4)  gives  a  brief 
development  of  the  equation  for  encounter  frequency  and 
for  the  case  of  ahead  seas,  the  incident  frequency,  is 


shown  to  be 


a)  =  3_ 

i  2V 


F¥, 


0)  -  1 
e 


This  relationship  was  plotted  by  Booth  and  is  repeated 
below  in  Figure  C-l.  For  the  XR-3  craft  operating  at  a 
steady-state  speed  of  30  knots,  V  =  50.67  ft/sec.,  the 
incident  wave  frequency  becomes 


a).  =  0.318 

i 


u 


1  +  (6.3), 


-i] 


An  approximation  to  this  functional  relationship  that  is 
useful  over  the  range 

1  <  a)  <20  rps 


2  =  -  /2  for  the  30  knot  speed  fo  the  XR-3  craft, 


IS  CO  .  =  (0 

i  e 
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* 


I 


♦ 


* 


Figure  C-l.  Wave  Frequency  vs.  Encounter 
Frequency 
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